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Abstract 



o 

$H Motivated from the theory of quantum error correcting codes, we investigate a combinato- 

Mh rial problem that involves a symmetric n-vertices colourable graph and a group of operations 

"^>( (colouring rules) on the graph: find the minimum sequence of operations that maps between 

^^ two given graph colourings. We provide an explicit algorithm for computing the solution of 

^\J our problem, which in turn is directly related to computing the distance (performance) of an 
underlying quantum error correcting code. Computing the distance of a quantum code is a 

^«>. highly non-trivial problem and our method may be of use in the construction of better codes. 

u 

r^' 1 Introduction 

-(— > 

a 

r^ With the recent advances in quantum information processing, protecting a physical system from 

I I environmental noise became an ultimate desideratum in the construction of a large scale quantum 

computer. Quantum error correcting codes, introduced about a decade ago [Il[2], play an important 

^ role in this respect and allow the protection of fragile quantum computation against undesired errors 

on via the use of redundancy. 

j^ An important class of quantum error correcting codes are the stabilizer codes [3j , which are the 

\^ quantum analog of classical linear codes [1]. Most known quantum error correcting codes belong 

.^ to the class of stabilizer codes, which in turn are equivalent |5j to the so-called graph codes. The 

^^ latter, as their name states, are codes built using an undirected symmetric graph on n vertices. For 

f^ more details about their properties see [6l[71|8]. 

. . The performance of quantum codes is characterized by the distance of the code, a positive 

. ^ integer that specifies how many errors the code tolerates (see e.g. [9j for a detailed discussion). The 

^ classical analog of the quantum distance is the well-known Hamming distance f3]. A vast amount of 

work has been dedicated to constructing quantum error correcting codes with the highest possible 
distance, see [10] for a list of up-to-date codes and bounds on the distance. 

Given a quantum error correcting code, determining its distance is highly non-trivial and in 
general requires numerical methods. Inspired by this, in the current article we investigate a combi- 
natorial problem involving graphs and operations on graphs, whose solution provides the distance 
of a corresponding graph code. Our solution is algorithmic and, although scales exponentiall}|jwith 
the number of vertices of the graph, can nevertheless be of great use in the construction of quantum 
error correcting codes. The solution we provide is elegantly described by linear systems over finite 



^Up to our current knowledge, determining the distance of a quantum code seems to be a computationally-hard 
problem. 



fields with the properties of the solution space providing the distance of the code. Although the 
space of solutions itself can be compactly described using a suitable basis, computing the distance 
reduces to searching among all possible solutions, and is this searching that makes the problem 
exponentially- hard. 

Even though the problem we investigate is motivated from the theory of quantum error cor- 
recting codes, we believe, given its combinatorial nature, that it may be of intrinsic interest to the 
discrete mathematics community, see e.g. [llj for a related work. 



2 Graph Labellings and Diagonal Distance 

Let Gn = (y, E) be a graph with n vertices where V is the vertex set of Gn and E is the edge 
set of Gn- A graph labelling of G„ is a function L : V ^ Z/2Z, which can be written as an 
element of (Z/2Z) by using the ordering of the vertices. For example, if we have a graph with 
V = {vi,V2,V3,V4,V5} and a graph labelling L = {(fi,0), (t;2, 1), (ws, 0), (f4, 0), (fs, 1)}, then we 
write the graph labelling simply as (0, 1,0,0, 1) G (Z/2Z) by requiring that the i^^ entry in the 

element of (Z/2Z) is the label of the i^^ vertex of the graph for 1 < i < 5. In order to simplify 
notation we will henceforth refer to the graph labelling L as the element of (Z/2Z") instead of as 
the function L : V ^ (Z/2Z) . In order to define the diagonal distance of a graph Gn we need to 
define two operations which indicate the possible ways to move between graph labellings. Defining 
these operations requires the notion of the adjacency matrix F of a graph Gn which is defined as 
the n X n matrix with Fjj = 1 if {vi,Vj) G E (there is an edge between Vi and Vj) and Fjj = if 
{vi,Vj) ^ E (there is no edge between Vi and Vj). Observe that (Z/2Z) is the set of all graph 
labellings of a graph with n vertices and so our operations our defined as follows. 

Definition 1. Define two operations applied to the graph labelling L of Gn for arbitrary Vi £ V by, 

• Zi : (Z/2Z)" -^ (Z/2Z)" defined by Li ^ Li + 1 mod 2. 

• Xi : (Z/2Z)'' -^ (Z/2Z)" defined by Lj ^ Lj + Tij mod 2, VI < j < n. 
where Li denotes the i*^ entry in L. 




Figure 1: An example of the two operations over (Z/2Z) applied to a graph with 7 vertices. 



We name an element in the group of operations at the i vertex of G„ by Oj G {Zi, Xi); that is, 
Oi is an arbitrary element of the group of operations generated by the operations Zi and Xi with 
composition of maps as the group operation. The characteristic r/-function is defined by, 



vio^ 



1 iiOi^I 
if Oi = I 



(1) 



where / : (^jll}^ — t- {%jll}j is the identity map. We are now prepared to define the diagonal 
distance of a graph. Intuitively, the diagonal distance of a graph G„ is the minimum number 
of composed Oj's that change the labelling L of G„ to a non-zero number of other labellings 
before returning to the original labelling L. Formally, summing the characteristic 77-function over 
a collection of Oj's counts the number of composed Oj's different from the identity. By ensuring 
that the composition (denoted by J|) of the Oj's equals the identity map /, we guarantee that the 
graph returns to its original labelling. Notice that the diagonal distance is invariant under different 
numberings of the vertices, as long as the notation remains consistent throughout the usage of this 
algorithm. 

Definition 2. The diagonal distance A(G„) of a graph Gn is defined as, 

A{Gn) = mm\Y,v(.Oi)>0 llO^ = I,0^ € {Zi,Xi),l < i < n\ (2) 



, i=l 



i=l 



3 Computing Diagonal Distance Over Z/2Z 

Given a graph G„, we construct the n x 2n matrix A by concatenating I„ and T as, 

A = [/„ I r] (3) 

where /„ is the n x n identity matrix and T is the adjacency matrix of G„. We think of /„ as the 
part of A which keeps track of Zi applied to Gn and F as the part of A which keeps track of Xi 
applied to G„. Let k"^ be an arbitrary nonzero 1 x 2n vector with entries in Z/2Z and write 

k^ = [ki-- -knl kn+l ■ ■ ■ k2n] (4) 

If Ak = (mod 2), that is, if k is in the nuUspace of A over Z/2Z, then we claim that k contains the 
information about the application of Oj's to G„ which changes the labelling L of G„ to a non-zero 
number of other labellings before returning to the original labelling L. In order to justify this claim 
we require the characteristic ^-function defined for 1 < i < n hj 

[0 if fej = fcj+n = 

By relating the characteristic x-function and the characteristic ry-function to each other we can 
relate elements Oi £ {Zi,Xi) to a vector k that is constructed to encode the information in Oj. More 
precisely, we consider an arbitrary set of group operations {Oj}i<i<„ and for each Oi = Zf^X^^ 
construct the vector k with the i entry given by ki = Zi {mod 2) and the {i + n) entry given by 
ki-\-n = Xi {mod 2), where 1 < i <n. 

Lemma 1. Let {Oi}^^i he an arbitrary collection of Oi G {Zi,Xi). Then for any Oi in this 
collection, there exists Zi,Xi G Z/2Z such that Oi = Z^'-X^' and 

V{0i) = x{ki) (6) 

where k is the 2n x 1 vector with ki = Zi and ki^n = Xi for 1 < i < n. 



Proof. Let Oi = Z^'-X^^, where Zi,Xi € Z. Note that since labelhngs are elements of Z/2Z, then 
Zf = I and Xf = I since these maps act by adding the same quantity twice. Thus Zl^X^^ = 
Z^^ mo -^x^ mo ^^^ g^ Q^ jg equivalent to ZfXf for some x, z G Z/2Z. Construct the vector k 
as stated in the lemma. If r]{Oi) = 1, then ZfXf ^ I and so either z ^ or x ^ 0, which implies 
that xiki) = 1- If viOi) = 0, then ZfXf = I and so z = x = 0, which implies that xiki) = 0. 
Therefore, r]{Oi) = x{h)- D 

We now prove another lemma which restricts the possible vectors k, and thus the possible set 
of group operations {Oi}i<j<n) which can satisfy the definition of diagonal distance. 

Lemma 2. Let {Oj}i<j<n be a set of group operations with each Oi = Z^^X^\ where Zi,Xi G Z/2Z. 
Then, 

n 

YlOi = I <^^ke ker(A)/2Z (7) 

1=1 
where k is the 2n x 1 vector with ki = Zi and fcj+n = Xi for 1 < i < n. 

Proof. Write a labelling L G (Z/2Z) of a graph G„ as an n x 1 vector and represent {Zi,Xi) by 
additive matrix operations as Zi{L) = L + Ci mod 2 and Xi{L) = L + coli(r) mod 2, where e^ is 
the n X 1 matrix with 1 in the i**^ entry and elsewhere and where colj(r) is the i**^ column of F. 
(=^) Then for an arbitrary Oi = Z^^X^^ we have that Oi can be represented by an action on the 
labelling L as Oj = L + ZjCj + Xicolj(r) mod 2. Then for any set of group operations {Oi}i<i<„, 

n n n 

Y[Oi = L + '^Ziei + '^XiColiiT) mod 2 = L + Ak mod 2, (8) 

i=l i=l j=l 

If nr=i Oi = I, then nr=i Oi = L mod 2 which yields 

L + Ak = L mod 2 (9) 

and therefore every set of group operations {Oi}i<i<„ has a representation as a vector in the kernel 
of A. 

(<^) The backward implication follows from rewriting the previous equation, since we now assume 
that k is some arbitrary vector in the kernel. Use Lemma 1 to construct a set of group operations 
corresponding to k. Since k G ker(A)/2Z by assumption, we get that 

L + Ak = L mod 2 

Then, as in the other direction, we have that nr=i Oi — L, which implies that HILi Oi = I- 
Therefore, 

n 

k G ker(A)/2Z ^^ JJ Oi = /. (10) 

i=l 

D 
This leads to the main theorem regarding the computation of diagonal distance. 
Theorem 1. The diagonal distance A(G„) of a graph Gn can he computed as 



A(G„)=mini^x(^*)>0 



. j=i 



k G ker(A)/2Z \ . (11) 



Proof. By Lemma [T] and Lemma [2] we have that, 

( n n 

A(G„) = mill \ Y^niOi) > J]Oi = /,0, G {Z,,X,), 1 < i < 



. i=l 



n 



4 = 1 



mm <^ ^ x(fci) > 



, i=l 



A; e ker(A)/2Z \ . 



(12) 

D 



An immediate consequence of this theorem is a deterministic algorithm for computing the 
diagonal distance A(G„) of any graph G„. To run this algorithm, one must first compute a basis 
for ker(A)/2Z. Then the entire kernel must be spanned to check each vector k for x(^)) ^iid the 
smallest x{k) (along with its corresponding vector, if extra information is desired) is stored. Once 
the algorithm terminates, the smallest of all the x{k)^ will have been stored and will be equal to 
the diagonal distance A(G„). 

4 Generalization of Diagonal Distance to Multigraphs over l^jpL 

We now generalize the diagonal distance of a graph over Z/2Z to the diagonal distance of a multi- 
graph (a graph that can have multiple edges between vertices) over Z/pZ, where p is prime. We 
define a graph labelling exactly as over Z/2Z and note that the set of all graph labellings is now 
(Z/pZ) . The operations Zi and Xi applied to a graph labelling is defined exactly as over l^lTL, 
except the adjacency matrix F of a multigraph G„ is defined as the n x n matrix with Tij = k 
if there are k occurrences of (v, ,Vj) £ E (there are k edges between vi and Vj ) and Tij = if 
{vi,Vj) ^ E (there is no edge between Vi and Vj). We can then define the group of operations 





Figure 2: An example of the two operations Zf and Xf applied to a multigraph with 7 vertices. 



Oi = {Zi,Xi), with our operations defined over Z/pZ in order to obtain the straightforward gener- 
alization of Theorem 1 to Z/pZ by writing mod p everywhere that mod 2 is written in the proofs 
of Lemma [T| Lemma [2] and Theorem [l] 



Theorem 2. The diagonal distance A(G„) of a graph Gn can he computed as, 



AiGn)=mm\Y,x{h)>0 



, i=l 



k G ker(A)/_pZ 



(13) 



5 Minimal Diagonal Distance of Codes Over l^j'pL 

Thus far we have considered the diagonal distance A(G„) over "LJIL of a graph G^ and, more 
generally, the diagonal distance A(Gji) over "LlpL of a multigraph Gn- In either of these cases, we 
built into the definition of diagonal distance that HI^i Oi = I where / : (Z/pZ) — t- {l^j-pl}^ is the 
identity map. Yet, we can extend this notion to requiring that HILi ^« ~ -^' where L : (Z/pZ) — t- 
(Z/pZ) is some specific map, e.g. the map sending the configuration to 1. 

Let {Ci, ...,Gk} be a set of graph labellings (i.e. codewords), where Ci = 0. We then want 
to determine the diagonal distance between any two graph labellings Gr and Gs in this set, the 
distance of the code, which can be defined as follows. 

Definition 3. For a multigraph Gn and a set of graph labellings {Ci, C2, ..., C^}, the diagonal 
distance Ars{Gn) between two graph labellings Gr and Gg is defined as, 

( n n '\ 

A^,(G„) = mini^r/(O,)>0 [] O, = L,„ O, G (Z^, X^), 1 < i < n I (14) 



.i=l 



i=l 



where Gn is assumed to be originally labelled by Gr and Lrs is the map sending the graph labelling 
Gr to the labelling Gs- 

Similarly to computing A{Gn), we can minimize Ars{Gn) by minimizing the summation of the 
characteristic ^'function of vector k' such that Ak' = Gr — Gg mod p. It is a straightforward 
generalization of Lemma [T] and Lemma [2] in order to obtain the analogue of Theorem [2] for distance 
of a code. 

Theorem 3. The generalized diagonal distance Ars{Gn) of a graph Gn can be computed as, 

A,.,(G„) = minJ^x(^i) >0 Ak' = Gr - Gs mod p I (15) 

where Gn is originally labelled by Gr- 

Note that the distance (5 of a code {Ci, C2, ..., Gk} is simply min{Ars(G„) | r,s £ {1, 2, -..k}} 

6 Computational Example 

We run the main algorithm on the 5-cycle graph in Z/2Z. It is well-known that the diagonal 
distance of the 5-cycle is 3, for example by the set of operations shown in Figure 3. This can be 
found directly using our algorithm. To compute the diagonal distance, we must first set up the 
matrix A = [/„ | P]. For the 5-cycle, this is a 5x 10 matrix which has many possibilities depending on 
how the vertices are numbered. However, diagonal distance is invariant under different numberings 
and so we will just pick one. 

100000100 1' 
10 10 10 
A= 0010001010 
10 10 1 
0000110010 



Figure 3: An example of three operations that illustrate the diagonal distance of the 5-cycle. 



Next we must compute the kernel of this matrix. Some computation gives that the kernel has the 
following basis vectors. 

fci = (0,1, 0,0, 1,1, 0,0,0, of 
^2 = (1,0, 1,0, 0,0, 1,0,0, of 
^3 = (0,1, 0,1, 0,0, 0,1,0, of 
/C4 = (0,0,1,0, 1,0,0,0, 1,0)"^ 
^5 = (1,0, 0,1, 0,0, 0,0,0, if 

To find the vector in the kernel with the least amount of I's in it, we must search through all the 
linear combinations of these basis vectors. It is easy enough to write a program to do this. For 
this example, it actually turns out that the linear combinations with the least amount of I's in 
them are exactly the basis vectors themselves. Notice that all of these vectors have 3 I's in them. 
Therefore, the diagonal distance of the 5-cycle graph is 3. 
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